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Abstract. This paper presents a consistent description of the formation and the subsequent evolution of gaseous 
planets, with special attention to short-period, low-mass hot-Neptune planets characteristic of p Ara-like systems. 
We show that core accretion including migration and disk evolution and subsequent evolution taking into account 
irradiation and evaporation provide a viable formation mechanism for this type of strongly irradiated light planets. 
At an orbital distance a ~ 0.1 AU, this revised core accretion model leads to the formation of planets with total 
masses ranging from ~ 14 M® (0.044 Mj) to ~ 400 Mq (1.25 Mj). The newly born planets have a dense core of 
~ 6 M 0 , independent of the total mass, and heavy element enrichments in the envelope, Mz.env/Alenv, varying 
from 10% to 80% from the largest to the smallest planets. We examine the dependence of the evolution of the born 
planet on the evaporation rate due to the incident XUV stellar flux. In order to reach a p Ara-like mass (~ 14 M 0 ) 
after ~ 1 Gyr, the initial planet mass must range from 166 M 0 (~ 0.52 Mj) to about 20 M 0 , for evaporation rates 
varying by 2 orders of magnitude, corresponding to 90% to 20% mass loss during evolution. The presence of a core 
and heavy elements in the envelope affects appreciably the structure and the evolution of the planet and yields 
~ 8 % — 9% difference in radius compared to coreless objects of solar composition for Saturn-mass planets. These 
combinations of evaporation rates and internal compositions translate into different detection probabilities, and 
thus different statistical distributions for hot-Neptunes and hot-Jupiters. These calculations provide an observable 
diagnostic, namely a mass-radius-age relationship to distinguish between the present core-accretion-evaporation 
model and the alternative colliding core scenario for the formation of hot-Neptunes. 
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1. Introduction 

Since the discovery of the first exoplanet by Mayor & 
Queloz (1995), planet hunters have been discovering plan¬ 
ets with smaller and smaller masses. With the discovery 
of objects in the range of 10-20 Mq, i.e. 0.03-0.06 Mj ^ 
(Butler et al. 2004; McArthur et al. 2004; Santos et al. 
2004), a new step in the quest for Earth-like planets has 
been taken. While representing an extraordinary achieve¬ 
ment from the observational side, theorists still struggle 
to understand the structure and the origin of these light 
giant planets. Are they essentially composed of ices and 
rocks, with possibly a thin atmosphere, like our ice giants 
Uranus and Neptune? Or do they originate from larger 
gaseous planets, with a large gaseous envelope and a rela¬ 
tively small central rocky core? The answer to these ques¬ 
tions requires an understanding of their formation process. 
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Current planet formation scenarios, based either on 
the core accretion model or on gravitational instability, 
can more or less explain the presence of relatively massive 
planets with masses > 100 M 0 (or even larger in the case 
of the gravitational instability scenario) at various orbital 
separations but they do not necessarily predict the forma¬ 
tion of a large number of lighter planets (Boss 2001; Ida 
& Lin 2005; Papaloizou & Nelson 2005). 

In the framework of the core accretion model (Pollack 
et al. 1996), the general expectation was to find preferen¬ 
tially planets either less massive or more massive than the 
newly discovered Neptune-like planets, for these latter lie 
within the domain of critical mass (~ 10 — 20 Mq) above 
which runaway accretion of gas begins. Similarly, Ida & 
Lin (2004) suggest a possible deficit of intermediate mass 
planets (~ 10 — 100 Mq, i.e. 0.03 — 0.3Mj) with orbital 
separation a < 3 AU. On the opposite, using N-body simu¬ 
lations of colliding cores in a protoplanetary disk, Brunini 
& Cionco (2005) find that Neptune-like planets close to 
their host-star can form easily as a by-product of plane- 
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tary formation. If their scenario is correct, these planets 
should be of ice-rock composition with only a thin atmo¬ 
sphere. They also predict that a large population of these 
’’hot cores” should be discovered in a near future. 

As an alternative to this colliding core scenario, and 
given the close orbital distance of the Neptune-like planets 
discovered up to now, one cannot exclude the possibility 
that they formed initially as larger giant planets wich have 
undergone atmospheric evaporation during their lifetime 
(Baraffe et al. 2005). Within this picture, formation and 
evolution are strongly correlated: a correct understanding 
of light planet properties thus requires a consistent de¬ 
scription of the planet formation and evolution in order 
to interpret present-day observations. 

This paper is a first attempt to derive such a consis¬ 
tent picture from the planet formation to its subsequent 
evolution. We apply the core accretion models (cf. §2) de¬ 
veloped recently by Alibert et al. (2005a) to the 14 M 0 
(0.044 Mj) planet (modulo sin i) orbiting around the G- 
star fj, Ara at an orbital distance a = 0.09 AU (Santos et 
al. 2004). We focus on ^ Ara-like planets because it is still 
reasonable to apply the Alibert et al. (2005a) formation 
model at this orbital distance, whereas it is no longer the 
case at the location of the two other Neptune-like plan¬ 
ets which are located at a < 0.05 AU from their parent 
star (Butler et al. 2004; McArthur et al. 2004). Indeed, for 
a <0.1 AU, the description of the inner part of the disk, 
including tidal and magnetic interactions with the star, is 
too crude to provide a reliable formation scenario. 

Adopting the main characteristics of the newly born 
planets as predicted by the Alibert et al. (2005a) forma¬ 
tion model (core mass, heavy element content), we follow 
the later evolution of these planets, according to Baraffe et 
al. (2004, 2005), taking into account irradiation and evap¬ 
oration effects due to the vicinity of the parent star. We 
examine the sensitivity of the results on the evaporation 
rate by exploring a range of different rates and present our 
results in §4. Predictions and uncertainties of our scenario 
are discussed in §5. 

2. Formation model 

The adopted formation model is an extension of the core 
accretion model for giant gaseous planets developed by 
Pollack et al. (1996). The present model includes the ef¬ 
fect of migration and of disc evolution on the planet for¬ 
mation. Details can be found in Alibert et al. (2005a). 
The model describes successfully the formation of Jupiter- 
like planets on timescales consistent with typical disk life¬ 
times (Alibert et al. 2004). When applied to the formation 
of Jupiter and Saturn, the model yields core masses and 
heavy element contents in the envelope in agreement with 
recent interior structure models for Jupiter and Saturn 
(Saumon & Guillot 2004). Moreover, the model predicts 
surface compositions in agreement with both in situ abun¬ 
dance measurements by the Galileo probe for Jupiter, and 
with Earth based measurements for Saturn (Alibert et al. 
2005b). 


We employ this formation model to determine the 
properties of an ensemble of planets which can form 
eventually at 0.1 AU of a solar type star, like the 
planet found around /i Ara. We performed a large 
number of simulations, similar to those described in 
Alibert et al. (2005a, see their §3.2), exploring differ¬ 
ent initial parameters (orbital separation, dust-to-gas 
ratio in the disc, photoevaporation rate, disc mass). 
The initial disk profile is assumed to be a power law 
with a gas surface density E = Eo(r/5.2A[/)“^/^, with 
So = 100,300,500,700 g/cm^. The surface density of 
planetesimals is Sp = /^j/gS, where the dust-to-gas ratio 
Id/g is equal to 1/30 for temperatures below 150 K, and 
1/120 otherwise. These values reflect the high metallicity 
of the central star in the /i Ara system. Different initial 
orbital separations of the planet are adopted: Oini = 
0.2,0.3,0.4,0.6,1,1.3,1.6,1.9,2.2,2.5,3,3.5,4,4.5,5, 6 , 7 
AU. The evolution of the disk gas surface density is calcu¬ 
lated according to the a prescription (Shakura & Sunyaev 
1973), with a = 2 ■ 10“^, including a photoevaporating 
term (see Alibert et al. 2005a). Different values for the 
total photoevaporation rate are assumed: lO“®M 0 /yr, 
2 X lO-^Me/yr, 4 x lO-^Me/yr, 8 x lO-^Me/yr and 
15 X lO“®M 0 /yr. Finally, we assume a time delay between 
the beginning of disk evolution and the planet formation 
process. We adopt different time delays: 0 Myr, 0.5 Myr, 
1 Myr, 1.5 Myr, 2 Myr, 3 Myr, 4 Myr, 6 Myr and 8 Myr. 
The simulations include migration with a migration rate 
for low mass planets (type I migration, see Ward 1997) 
based on the results of Tanaka et al. (2002), reduced 
by a constant factor equal to 0.01 (see Alibert et al. 
2005a). As shown in Alibert et al. (2005b) for the cases 
of Jupiter and Saturn, a variation of the type I migration 
rate does not affect the final structure of the planet. Each 
simulation produces a planet, eventually falling onto the 
central star, and we focus on the properties of planets 
ending their formation around 0.1 AU, whatever their 
final mass. The results are summarised in Figure 1 which 
displays the core mass Mcore and the mass fraction of 
heavy elements in the envelope as 

a function of the total mass of the planet. Mz^env is the 
amount of mass deposited by the infalling planetesimals 
in the envelope whereas Mcore is the mass of solids 
that directly reach the central solid core of the planet. 
The mass of the core depends on some poorly known 
quantities, such as the characteristic size of the infalling 
planetesimals (we assume here a size of 100 km), and 
their internal properties (e.g. tensil strength, for which 
we use here values relevant for water ice). Note, however, 
that the total content of heavy elements is much less 
sensitive to these quantities. As shown in §4.2, the 
effect of evaporation depends mainly on the mean heavy 
element enrichment of the planet, and not on their precise 
distribution in the planet. Finally, we note that we do 
not take into account in the formation calculations the 
modifications of opacity and equation of state due to 
changes of the chemical composition in the envelope (see 
discussion in §5.2). 
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Fig. 1. Properties of planets forming at 0.1 AU from their 
sun, based on a set of simulations using different initial 
conditions (see text, §2). Lower panel: Core mass as a 
function of total mass. Upper panel: Metallicity in the en¬ 
velope .^env) defined as the ratio of heavy element mass in 
the envelope Mz,env to the total envelope mass Menv, as 
a function of total mass. 

The simulations predict that planets reaching a total 
mass in the range 14-400 M0 (~ 0.04—1.25 Mj) form with 
essentially the same core mass Mcore — SMq (cf. Fig. 1). 
The heavy element content in the envelope is found to 
increase substantially with decreasing total mass. The en¬ 
velope metallicity mass fraction Zenv varies from 80% for 
a 14 M® planet to 10% for masses > 100 M®. Note that 
our calculations do not take into account the probability 
of occurrence of the initial parameters leading to the for¬ 
mation of these planets. Therefore, the number of light 
planets that can actually form by this extended core ac¬ 
cretion mechanism and its fraction compared to the num¬ 
ber of more massive planets cannot be inferred from these 
calculations. 

3. Evolutionary models 

3.1. Irradiation and evaporation effects 

Using the initial conditions provided by the planet for¬ 
mation model described in §2, core mass and heavy el¬ 


ement enrichment for a total planet mass, we calculate 
the subsequent evolutionary sequences for different val¬ 
ues of these initial planetary masses. Evolutionary mod¬ 
els take into account irradiation effects from the parent 
star on the planet atmospheric and internal structure, 
as described in Baraffe et al. (2003) and Chabrier et al. 
(2004). The incident stellar flux received by the planet 
atmosphere is determined from the orbital and stellar pa¬ 
rameters characteristic of the /i-Ara system (a = 0.09 AU, 
TO* = 1 M®, ii* ~ 1 i?Q, Tefi* ~ 5800 K, cf. Santos et al. 
2004). Evaporation effects are also included, as described 
in Baraffe et al. (2004), using, as our fiducial model, the 
energy-limited escape model derived by Lammer et al. 
(2003, hereafter L03). This model describes a hydrody¬ 
namic mass loss process due to the high energetic XUV 
flux of the parent star. We also consider smaller evapo¬ 
ration rates. Indeed, L03 assume that the planet under¬ 
goes maximal energy-limited evaporation, thus providing 
an upper limit for escape rates. Their predicted escape rate 
for HD 209548b is ~ 100 times larger than the lower limit 
measured by Vidal-Madjar et al. (2003) for this planet. In 
a recent study, Yelle (2004) included a better treatment 
of atmospheric chemistry, a crucial ingredient for handling 
cooling processes, and found escape rates 20 or 100 times 
smaller than the ones estimated by L03. Taking into ac¬ 
count a two-dimensional, hydrodynamical energy deposi¬ 
tion calculation instead of the single-layer heating model 
used in L03, Tian et al. (2005) find escape rates for HD 
209458b at least ^ 16 times smaller than L03. Finally, in a 
recent study, Jaritz et al. (2005) analyze Roche-lobe effects 
on expanded upper atmospheres of close-in giant planets. 
They conclude that, in some cases, planets can undergo 
geometrical blow-off, rather than classical hydrodynamic 
blow-ofU, because the critical level where blow-off occurs 
cannot be reached before the exobase level reaches the 
Roche-Lobe. Jaritz et al. (2005) show that the transit 
planet OGLE-TR 26b is in such a configuration and And 
an escape rate 25 times lower than the maximal energy- 
limited mass loss calculated by L03. Given the large un¬ 
certainties in the escape rate illustrated by these different 
studies, we performed evolutionary calculations with 1, 
1/20 and 1/100 times the escape rate based on the model 
of L03, corresponding to values ranging from 5x10^^ g/s 
(2.6 10“® M®/yr) to 7x10® g/s (4 10“^^ M®/yr) for plan¬ 
ets older than 1 Gyr. 

3.2. Models with a rocky core 

All our previous evolutionary calculations were done for 
coreless planets (Baraffe et al. 2003, 2004, 2005, Chabrier 
et al. 2004). In order for the evolution to be consis- 


^ Classical dynamical blow-off occurs above a critical temper¬ 
ature for which the thermal velocity of atoms or molecules in 
the planetary upper atmosphere exceeds their escape velocity. 
In contrast, the geometrical blow-off, as deflned by Lecavelier 
des Etangs et al. (2004), occurs when the upper atmosphere 
reaches the planet Roche Lobe. 
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tent with our formation model, we have now included 
a central rocky core in our structure and evolution cal¬ 
culations. Instead of using a constant core density, as 
done in some exo-planet models (e.g Bodenheimer et al. 
2003), we have implemented the ANEOS equation of state 
(Thompson & Lauson 1972). This equation of state (EOS) 
describes the thermodynamic properties of different mate¬ 
rial of planetary interest (ice, dunite, iron), derived from 
the Helmholtz free energy. Even though this EOS has 
physical limitations^, it gives a thermodynamic descrip¬ 
tion of phase transitions and thermodynamic quantities 
relevant for the evolution of planets as well as Rosseland 
opacities and electron conductivities. In the evolutionary 
code, we integrate the structure equations from the cen¬ 
ter to the surface and at the core/envelope boundary, the 
EOS switches from ANEOS to the Saumon-Chabrier EOS 
(Saumon, Chabrier & VanHorn 1995). At this boundary, a 
density jump is present due to the change in composition 
but continuity in pressure and temperature is enforced. 
In order to check our implementation of ANEOS, we have 
compared the mass-radius relationships obtained for water 
ice and olivine (or dunite, Mg 2 Si 04 ) planets to the results 
displayed in Guillot (2005) and Guillot et al. (1996) based 
on the same EOS. We find an excellent agreement with 
these studies. 

As suggested by the formation model (cf. Eig. 1), we 
adopt in the following a core mass of 6 M 0 independent 
of the initial planet mass. We assume that the core is 
made of dunite, as representative of rock. This material 
yields typical mean densities in the core ~ 6-7 g cm“^. 
We also performed comparative calculations with water 
ice cores, corresponding to a lower mean density ~ 3 g 
cm“^. Adopting ice or dunite for the core composition 
slightly changes the mass-radius relationship for planets of 
identical core and total mass. However, this change has no 
effect on the main conclusion of our paper and we therefore 
adopt dunite cores in the following. 

3.3. Heavy element enrichment in the envelope 

As shown in Eig. 1, heavy element enrichment in the en¬ 
velope due to accretion of planetesimals during the for¬ 
mation process can be significant, up to 80% in mass for 
the lightest planet formed (see Alibert et al. 2005a for 
details on the calculation of the interaction between in¬ 
falling planetesimals and the envelope, and the resulting 
heavy element enrichment). Under the conditions of tem¬ 
perature and pressure characteristic of planet interiors, 
the contribution to the EOS due to the interactions be¬ 
tween particles is of the order of or even dominate the 
kinetic (perfect gas) contribution. However, an accurate 
treatment of the EOS of such an interacting multi-species 
fluid is elusive and approximations are necessary. Since 
our goal is to explore the effect of a heavy element con¬ 
tamination on the planet structure and evolution, we treat 

® e.g electrons are treated classically, i.e no degeneracy; see 
also the discussion in §3.3.4 of Guillot et al. 2004 


this contribution to the EOS as follows. The presence of 
heavy elements with a mass fraction Z in the envelope is 
mimicked by an equivalent helium fraction Uequiv = Y + Z 
in the H/He Saumon-Ghabrier EOS, where Y is the gen¬ 
uine helium mass fraction in the envelope. This procedure 
has been discussed in Ghabrier et al. (1992) for models 
of Jupiter and Saturn and in Chabrier & Baraffe (1997), 
who demonstrate its validity for solar metallicity i.e. for 
Z <C 1. It is clear that for metal enrichment as high as 
obtained by the formation model, such an approximation 
becomes extremely crude. Since the heavy element enrich¬ 
ment will dominantly affect the structure of the planet 
(the thermal contribution from the solid core is small) and 
since within first order the different contributions to the 
EOS scale with the density, this approximation, however, 
is good enough to give a semi-quantitative description of 
the effect of an increase of Z on the thermodynamic prop¬ 
erties of the envelope material, and thus on the mechani¬ 
cal properties (mass-radius relationship and, hence, mean 
density) of our planets. 

4. Results 

We followed the evolution of planets of different initial 
masses, mi, in the range ~ 18-175 Mq (0.055 to 0.55 Mj). 
The initial interior structures Mcore and Zenv for the evo¬ 
lutionary models are the ones obtained from the formation 
model. We adopt an initial radius which is arbitrarily large 
but lies well within the Roche lobe radius of the planet. 
Work is in progress to better determine the final radius 
of the new born planet once accretion has terminated in 
order to use it as a more accurate initial condition for the 
evolutionary models. It should be stressed, however, that 
no planet formation model can provide an accurate de¬ 
termination of the final planet radius, at the end of the 
gas accretion phase. Indeed, the calculation of this final 
radius is a daunting problem, which depends crucially on 
the accretion process, not expected to be spherical, and 
on the ability of the accreting planet to radiate away its 
thermal energy, which requires a correct description of the 
(accretion) radiative shock. According to the results of the 
formation model, we assume that all planets have a rocky 
(dunite) core mass of 6 Mq and a heavy element enrich¬ 
ment corresponding to their mass (see Fig. 1). Evaporation 
rates are varied between 1 and 1/100 times the rate given 
by L03. We focus on sequences which, for a fixed evapora¬ 
tion rate, reach a mass close to a Neptune-mass (17 M 0 
or 0.053 Mj) within a few billion years, the inferred age 
for p-Ara (Santos et al. 2004; Baraffe et al. 2005). Tabled 
summarizes the different evolutionary sequences presented 
below, with their inputs and parameters. 

4.1. Evolution with maximal escape rate 

In the case of the maximal evaporation rate (rate from 
L03,) the initial masses must be greater than about 0.5 Mj 
(160 M 0 ) to produce a p-Ara like Neptune-mass planet af¬ 
ter a few billion years of evolution, and > 0.3 Mj, about 
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Table 1. Properties of evolutionary sequences with different initial masses, escape rate prescriptions, core masses and 
envelope metallicities, which can reach a mass close to the Neptune mass in more than 1 Gyr. 


Initial mass 

escape rate 

Mcore/ M® 

Zenv 

log t (yr) 

M/M® (Mj) 

R/Rj 

M (M®/yr) 

166 M® 

LOS 

0 

0.02 

7.0 

166. (0.52) 

1.44 

3.90 10~'^ 





8.0 

1S7. (0.43) 

1.28 

3.32 10“'^ 





9.0 

32. (0.10) 

1.25 

1.18 10“^ 





9 . 0 s 

20. (0.06) 

1.40 

2.45 10“^ 

166 M® 

LOS 

6 

0.02 

7.0 

166. (0.52) 

1.39 

3.48 10~" 





8.0 

140. (0.44) 

1.23 

2.85 10-'^ 





9.0 

70. (0.22) 

1.08 

3.57 10“® 





9.44 

21. (0.07) 

0.94 

2.92 10“® 

166 M® 

LOS 

6 

0.1 

7.0 

166. (0.52) 

1.33 

3.09 10“' 





8.0 

143. (0.45) 

1.18 

2.47 10“'^ 





9.0 

86. (0.27) 

1.05 

2.69 10“® 





9.80 

18. (0.06) 

0.83 

1.07 10“® 

50 M® 

LOS/IO 

6 

0.2 

7.0 

50. (0.15) 

1.52 

1.5 10"' 





8.0 

41. (0.13) 

1.19 

8.83 10“® 





9.0 

23. (0.07) 

0.93 

6.84 10“® 





9 . 2 s 

19. (0.06) 

0.87 

3.99 10“® 

33 Mq 

L0S/20 

6 

0.4 

7.0 

33. (0.10) 

1.31 

7.30 10"® 





8.0 

28. (0.09) 

1.08 

4.79 10“® 





9.0 

20. (0.06) 

0.82 

2.72 10“® 





10.0 

15. (0.05) 

0.64 

2.02 10"“ 

18 Mq 

LOS/lOO 

6 

0.4 

7.0 

18. (0.06) 

1.35 

3.01 10~® 





8.0 

16. (0.05) 

1.03 

1.43 10“® 





9.0 

14. (0.04) 

0.76 

6.02 10"“ 





10.0 

13. (0.04) 

0.61 

4.02 10"“ 


a Saturn-mass, to reach this limit after about 100 Myr. 
Figure|2shows evolutionary sequences of planets with ini¬ 
tial mass mi=0.3 Mj and 0.52 Mj. The results for mi=0.52 
Mj (166 M 0 ) illustrate the significant effect of the pres¬ 
ence of a core and of the envelope metal enrichment on the 
planet evolution. Both effects yield an increased average 
metallicity and thus an increase of the mean density of 
the planet. In order to highlight the effect of the presence 
of a central core and heavy element enrichment on the 
structure of a planet. Table |21 compares, for constant mass 
evolution, i.e. without evaporation, the radius of a 0.3 Mj, 
i.e. Saturn-mass, coreless gaseous planet of solar compo¬ 
sition and the one of a Mcore = 6 M 0 and Zenv = 0.1 
planet of same mass. The effect of such an internal com¬ 
position translates into a ~ 8 % — 9% difference in the 
radius for this planet mass, whereas for Jupiter-mass or 
larger planets, the effect is < 5%. Inferring the presence of 
significant heavy-element enrichment in the interior of ex¬ 
oplanets should thus be in reach for transit planets smaller 
than about a Saturn-mass. Figure 13 displays the interior 
pressure and density profiles of a TOi=0.52 Mj planet, with 
or without core and/or heavy element enriched envelope. 

Since the evaporation rates is inversely proportional to 
the mean planet density, m cx {Gpp)~^ (see Baraffe et al. 
2004, Eq. (1)), they are affected by the presence of a core 
and an increase of Zenv For the sequence with mi=0.52 
Mj and with the structure predicted by the formation 
model (Mcore =6 M 0 and Zenv=0.1, solid line in Fig.EJ, 
the sequence reaches masses < 20 M 0 (0.06 Mj) within 


Table 2. Effect of heavy element content on the radius 
of a 0.3 Mj (~ 100 M 0 ) at different ages. Evolution 
proceeds at constant mass, i.e M = 0. Case A corre¬ 
sponds to Mcore = 0, Zenv=0.02; Case B corresponds to 
Mcore = 6 M 0 , Zenv=0.1. Irradiation effects are included 
in both cases. 


age 

R/Rj 

R/Aj 

(Gyr) 

(case A) 

(case B) 

0.5 

1.19 

1.09 

1 

1.14 

1.05 

5 

1.06 

0.98 


6 Gyr with an escape rate ~ 2x 10^^ g/s ( 10“® M 0 /yr) at 
this age (see table Without a core and heavy element 
enrichment, this mass domain is reached within ^ 1 Gyr. 
Fig.Elalso reveals that when adopting the maximal evapo¬ 
ration rate, planets with initial masses ^ 100 M 0 (0.3 Mj) 
reach the runaway phase described in Baraffe et al. (2004) 
and evaporate entirely within ~ 100 Myr regardless of the 
presence or not of a core and heavy element enrichment. 
We must also recall that the evaporation rates depend on 
the high energy flux from the star, which decreases ex¬ 
ponentially with time. This decrease, combined with the 
rather low average density of planets at early ages, yields 
escape rates at the beginning of the evolution which are 
^ 30 times higher than those obtained after 1 Gyr (see 
Baraffe et al. 2004). 
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7 8 9 10 

log t 

Fig. 2. Evolution of irradiated and evaporating planets at 
an orbital distance a=0.09 AU from a G-type star. The 
escape rate of LOS is used. Sequences with initial mass ~ 
166 M0 (0.52 Mj) are shown for different cases: Mcore=0 
and metallicity in the envelope ^env=^o (dot); Mcore=6 
M0 and Zenv=^0 (dash); Mcore=6 M0 and Zenv=0.1 
(solid). The (red) dash-dotted curve corresponds to an 
initial mass ~ 100 M0 (0.3 Mj), with Mcore=6 M0 and 
^env=0.15, as suggested by Fig. 1. 

4.2. Evolution with lower evaporation rates 

We now consider the effect of different evaporation rates, 
below the LOS value, on the evolution of the planets. Since 
the present formation model (§2) can produce planets with 
masses smaller than ~ 100 M0 (0.3 Mj), we consider var¬ 
ious objects below this limit. Figure 0 portrays the evolu¬ 
tion of planets with initial masses of 50 M0 (~ 0.15 Mj) 
and 33 M0 (~ 0.1 Mj), respectively, for evaporation rates 
corresponding to 1/10 and 1/20 times the LOS value. For 
the corresponding planet densities, these rates correspond 
to M <10^^ g/s (6.8 ■ 10“® M0/yr) and <5 x 10^^ g/s 
(2.7-10“® M0/yr) after 1 Gyr, respectively (see tabled]). 
As mentioned previously, such 10 to 20 smaller rates than 
the LOS maximum escape value are suggested by various 
calculations (Lecavelier des Etangs et al. 2004; Yelle 2004) 
and in particular by the recent multi-layer hydrodynami- 
cal calculations of Tian et al. (2005). The planets have a 6 
M0 core and a heavy element enrichment in the envelope 
of Zenv=20% and Zenv=40% in mass fraction for initial 
masses 50 M0 and 33 M0 respectively (Fig. 1). As seen 



2 I I I I I I I I I I I I I 1 I I I I I I i 

0 0.2 0.4 0.6 0.8 1 


r/Rp 

Fig. 3. Pressure (upper panel) and density (lower panel) 
radial profiles (in units of total radius) for a planet with 
mass ~ 166 M0 (0.52 Mj). The structures displayed cor¬ 
respond to an age of 1.5 Gyr. Different cases are shown: 
Mcore=0 and metallicity in the envelope Zenv=^0 (dot); 
Mcore=6 M 0 and ^env=^0 (dash); Mcore=6 M 0 and 

^env = 0.1 (solid). 


in Fig.^ for such rates, the planets reach a Neptune-mass 
within a few hundred of million years to a few gigayears, 
meaning that about 2/3 to 1/2 of their initial mass has 
been lost under the influence of the parent star high en¬ 
ergy flux. 

If we adopt an escape rate corresponding to 1% of the 
L03 value, i.e. a reduction by a factor 100, in order for 
the planet to reach a Neptune-mass within ^ 1 Gyr, the 
initial mass must be rrii ~ 18M0 (0.056 Mj), as illus¬ 
trated in Fig. El This corresponds to a ~ 20% mass loss of 
the planet envelope. For this initial mass, our formation 
models predict an envelope enrichment .^env — 0.4. 

This abundance of heavy elements strongly depends on 
the treatment of planetesimal accretion (size of the plan- 
etesimals, mechanical destruction in the envelope, convec¬ 
tion in the envelope, see §5.2). To test how sensitive our 
evolutionary models are to these issues, we performed a 
comparison calculation in which we assumed that all the 
accreted planetesimals sink through the envelope without 
mass loss and fall onto the core. In such a model, a planet 
with an initial mass of 18 M0 is found to have a core mass 
of ^ 11 M0 surrounded by an envelope essentially com¬ 
posed of H and He. In that case, the evolutionary proper¬ 
ties are very similar to the ones obtained with our nominal 
model (characterized by Mcore=6 M0 and Zenve=0.4) for 
L03 evaporation rates divided by 100. The radius found in 
the former case is ~ 0.55 i?j at 10 Gyr, instead of 0.6 i?j 
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Fig. 4. Same as Fig. El for planets with initial mass 50 
M0 (0.15 Mj) and 33 M© (0.1 Mj) and escape rates 1/10 
(dash) and 1/20 (solid) times the rates of L03. For both 
cases, Mcore=6 M0. For the 50 M0 planet (upper dashed 
curve in both panels), Zenv=0.2 and for the 33 M0 planet, 
Zenv=0.4 (cf. Fig. 1). 

for the nominal model (solid curve in Fig.^J. The different 
evolutionary sequences are summarized in Tabled 

5. Predictions and uncertainties 

5.1. Predictions 

Our results predict that if the formation and evolution of 
a hot-Neptune, such as the planet around /i-Ara, proceeds 
as described by one of the sequences displayed in Figs.EJ- 
El its radius should be greater than ~ 0.6 i?j, regardless 
of the escape rate adopted. Baraffe et al. (2005) reached 
similar conclusions for coreless evolutionary sequences and 
maximal escape rate. The present work confirms this pre¬ 
diction with more realistic planetary structures, consistent 
with the core accretion model scenario. It is important to 
understand the reason for such a large radius and thus 
to disentangle the various effects, evaporation, irradiation 
and heavy element enrichment. This is illustrated in Table 
Elwhich displays a few values of the radius of Neptune mass 
planets evolved at constant mass with different Zenv and 
with or without irradiation effects included. We first recall 
that evaporation effects do not affect the evolution, and 
thus the mass-radius relationship of the planet as long as 
this latter is not in the evaporation runaway phase (see 
Baraffe et al. 2004). The values of radii given in Table El 
are thus also characteristic of radii obtained for evolution¬ 
ary sequences including evaporation for same mass and 
age (if they are not in the runaway phase). As seen in 



7 8 9 10 

log t 


Fig. 5. Same as Fig. El for planets with initial masses 18 
M 0 (0.056 Mj), Mcore=6 M 0 and Zenv=0.4 (see Fig. 1) 
but different escape rates: 1/20 (dash) and 1/100 (solid) 
times the rates of L03. 


the table (case 1 vs case 2), an increase of heavy element 
from Zenv=0.1 to Zenv=0.4 yields a ~ 15% decrease of 
the radius of an irradiated Neptune mass planet at 5 Gyr 
while, for a given metal enrichement Zenv=0-4, irradia¬ 
tion effects increase the radius by ^ 20-30% at a given 
age. Table 0 also shows that radii as large as 0.8 i?j after 
a few Gyr can be reached for an envelope metal fraction 
Zenv=0.1 in the irradiated case. This envelope metallicity 
is predicted by our formation model if the initial planet 
mass is > 150 M0 (~ 0.6 Mj). In order for such a high 
progenitor mass to reach a Neptune mass planet within 
a few gigayears, high escape rates are required (see §4.1 
and Fig. El This analysis shows that the aforementioned 
~ 0.6 i?j lower limit for the radius of hot Neptunes stems 
essentially from the effect of irradiation on a planet which 
retains a substantial gaseous (H, He) envelope. 

In contrast, if /r-Ara-like planets are hot cores with 
small gaseous envelope, as suggested by Brunini & Cionco 
(2005), their radius should be significantly smaller, closer 
to the Neptune planet radius (~ 0.35 Rj). Radius deter¬ 
minations of hot-Neptunes could thus distinguish between 
different formation scenarios. They may also provide in¬ 
formation about the efficiency of the evaporation process, 
since measurements of radii larger than ~ 0.8 Rj, for plan¬ 
ets older than a few Gyr, would indicate high initial pro¬ 
genitor mass and thus high escape rates, as above men¬ 
tioned. Evaporation rates also affect the time spent in the 
mass range 10-20 M0 (0.03-0.06 Mj). Evolution of plan¬ 
ets at the maximal evaporation rate proceeds rapidly and 
the corresponding detection probability while the planet 
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Table 3. Effect of heavy element content and of irradia¬ 
tion on the radius of a Neptune mass planet 17 M 0 ) at 
different ages. Evolution proceeds at constant mass, i.e M 
= 0. All models have Mcore = 6 M 0 . Case 1 corresponds 
to Zenv=0.1, with irradiation effects (from a Sun at 0.09 
AU); Case 2 is similar to case 1 with Zenv=0.4. Case 3 is 
with Zenv=0.4, without irradiation effects. 


age 

R/Rj 

R/Rj 

R/Rj 

(Gyr) 

(case 1) 

(case 2) 

(case 3) 

1 

0.93 

0.80 

0.62 

5 

0.83 

0.70 

0.58 


lies within this mass range is much smaller than for sig¬ 
nificantly reduced rates. 

The internal composition also strongly affects the evo¬ 
lution of the planet. The sequence with initial mass 0.52 
Mj (166 M 0 ), with Mcore=6 M 0 and Zenv=0.1 (solid line 
in Fig. 121), spends 1.1 Gyr with a mass < 30 M 0 (0.094 
Mj). For the same object without heavy element enrich¬ 
ment in the envelope (dashed line in Fig.[21), this time is 
reduced by a factor ^ 4. Finally, for the same planet with¬ 
out a core and heavy element enrichment (dotted line in 
Fig. 121), the time is reduced by a factor ^ 15. A planet 
with initial mass 1 Mj, with the same Mcore =6 M 0 and 
■^env=0.1, looses Only 15% of its total mass after 10 Gyr, 
using the rates of L03 at an orbital distance a = 0.09 
AU. Jupiter-like planets thus remain almost unaffected by 
evaporation at such orbital distances. Assuming the same 
formation probability for a 0.5 Mj and a 1 Mj planet, our 
model predicts for the high LOS evaporation rate a signif¬ 
icantly lower detection probability for hot-Neptunes than 
for hot-Jupiters. Roughly 10 times smaller, comparing the 
1.1 Gyr lifetime of a planet with mass < 30 M 0 to the 10 
Gyr lifetime of a Jupiter-mass planet. Note that without 
heavy element enrichment (in the core and in the enve¬ 
lope), this probability would even be signihcantly smaller. 

5.2. Uncertainties and model improvements 

In the following, we describe a number of assumptions 
used in the core accretion model and the evolutionary cal¬ 
culations, which need to be improved in the future. 

A main uncertainty in the formation model stems from 
the treatment of the accreted matter and its contribution 
to the gravitational potential of the planet. As done in 
Pollack et al. (1996), the mass of accreted planetesimals 
is systematically added to the core mass in the gravita¬ 
tional potential calculation, whether they are destroyed 
in the envelope or fall onto the core (this is the so called 
’’sinking approximation”, see Pollack et al. 1996, Alibert 
et al. 2005a). Work is underway to include a more consis¬ 
tent treatment of the gravitational potential, taking into 
account the distribution of the accreted matter within 
the envelope. Note, however, that planetesimals are de¬ 


stroyed very deep in the envelope"^, in regions close to the 
core/envelope interface. 

Goncerning evolutionary models, we have already dis¬ 
cussed the assumption used for the treatment of heavy 
element enrichment in the envelope (see §3.3). A more 
consistent treatment of such enrichment may affect the 
envelope structure and, hence, the mean density of the 
planet and the evaporation rates. It is difficult to esti¬ 
mate the differences between an envelope structure based 
on the assumption of an effective helium fraction, as used 
here, and the one based on a consistent EOS. For values 
up to .Z w 10 % of heavy element enrichment in the en¬ 
velope, however, we anticipate that a proper treatment 
of the equation of state should not affect drastically the 
present results. 

Another important assumption stems from the irradi¬ 
ated atmosphere models used as outer boundary condi¬ 
tion for the interior structure. These models (Barman et 
al. 2001, 2005) include irradiation effects with up-to-date 
molecular and grain opacities, but have been computed 
for solar composition. Taking into account the increase 
of metallicity in atmospheric calculations up to levels in¬ 
ferred by the present formation models is a complex task 
and a long term project (see Chabrier et al. 2006 for pre¬ 
liminary results). We thus use here solar composition at¬ 
mospheric models. Note however that for close-in planets, 
irradiation effects from the parent star yield planetary at¬ 
mospheric structures which are radiative even at deep lev¬ 
els (Barman et al. 2001; Baraffe et al. 2003). Therefore, 
when the proto-planet accretes, the bulk of planetesimals 
either settle down to the core or are destroyed in the deep 
envelope. Only the smallest planetesimals may pollute the 
atmosphere. If convection in the interior is efficient, heavy 
elements will be mixed homogeneously, but should not 
pollute the outermost irradiated, radiative surface layers. 
This is no longer true, of course, for non-irradiated plan¬ 
ets, as measurements of the Jupiter atmosphere by the 
Galileo probe indicate abundances of CH 4 , H 2 S, NH 3 three 
times solar. For hot-Neptunes and hot-Jupiters, we may 
thus expect the heavy element enrichment in the atmo¬ 
sphere to be less dramatic than in the convective envelope. 
This point is also important for the treatment of evapo¬ 
ration, since the model of L03 applies to hydrogen domi¬ 
nated upper atmosphere. The effect of heavy element con¬ 
tamination on evaporation rates still needs to be explored. 
Intuitively, one expects that an enrichement in heavy ele¬ 
ments yields a less efficient evaporation process. However, 
the effect is certainly more complex, since a change of the 
chemical composition will affect the whole chemistry and 
cooling properties of the upper atmosphere. The effects of 
heavy element enrichment on atmosphere structures and 
evaporation process are thus important problems to study 
in the future. 


^ deeper than in former calculations by Pollack et al. 1996, 
see Alibert et al. 2005a 
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6. Conclusion 

The present study presents first consistent calculations be¬ 
tween planet formation, within an improved version of the 
core accretion model, and subsequent evolution, including 
irradiation and evaporation. These calculations can ex¬ 
plain, among others, the existence of hot-Neptune planets, 
like the one recently discovered around ^ Ara. The forma¬ 
tion model constructed by Alibert et al. (2005a) allows the 
existence of planets with initial mass ranging from ~ 14 
M 0 (0.044 Mj) to ~ 400 M 0 (1.25 Mj). The planets end 
up having a dense core of ~ 6 M 0 , independently of the 
total mass, and heavy element enrichment in the envelope 
varying from 10 % to 80% from the largest to the smallest 
mass. With maximal evaporation rates, as predicted by 
the models of LOS, the progenitor of gigayear old fj, Ara- 
like planets has an initial mass ~ 0.5 Mj. For evaporation 
rates a factor 10 to 20 smaller, the progenitor has a mass 
in the range ~ 0.1-0.15 Mj (~30-50 M 0 ), meaning the 
planet has lost 1/2 to 2/3 of its initial envelope. For plan¬ 
ets already forming in the Neptune-mass range, < 20 M 0 , 
the LOS rate must be reduced by a factor 1/100 for the 
planets to survive to evaporation at 0.09 AU from their 
Sun. In this latter case, the evaporation rates vary from 
~ 10^^ g/s (6 10“^° M 0 /yr), at an age of 1 Gyr, to ~ 7 
10® g/s (4 10“^^ M 0 /yr) at 10 Gyr. Our current knowl¬ 
edge of evaporation processes in exoplanet atmospheres is 
still too embryonic to favor or exclude any of these values, 
although recent hydrodynamical simulations favor the in¬ 
termediate value (Tian et al. 2005). The range of escape 
rates explored in the present calculations, however, should 
bracket the ’’real” solution. 

Our calculations provide an observable diagnostic for 
such a core-accretion-irradiation scenario for the forma¬ 
tion and evolution of hot-Neptunes, namely radii > 0.6 
Rj for Neptune-mass planets. Also, the scenario with max¬ 
imal evaporation rate may provide a statistical signature 
in terms of the number of hot-Neptunes, which is expected 
to be much smaller than the number of hot-Jupiter-like 
planets. In the absence of detailed statistical analysis, a 
crude comparison of evolutionary timescales suggests 10 
times less hot-Neptunes than hot-Jupiters at a ~ 0.09 AU, 
for the LOS maximal escape rate. The detection of a larger 
fraction of hot-Neptunes at such orbital separation would 
thus favor scenarios with lower evaporation rates. 

A complete understanding of the formation and evolu¬ 
tion of light, Neptune-mass planets is still far from reach. 
The present attempt to derive a consistent picture from 
the planet genesis to its present day conditions provides 
one possible alternative path for the birth and fate of these 
objects. It also points to several issues to be addressed for 
a better modelling of planet formation and interior and 
atmospheric structures. Further improvement requires in 
particular (i) a better treatment of heavy element enrich¬ 
ment, both in the envelope and in the atmosphere and (ii) 
a better understanding of evaporation processes. 

In parallel with such theoretical developments, obser¬ 
vations of planets down to a few Earth-masses and lo¬ 


cated at orbital separations large enough (a > 1 AU) to 
remain unaffected by evaporation would provide strong 
constraints on the formation scenarios of light planets. 
This is another appealing challenge for planet hunters. 
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